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PARAMETER  LIST 


x,y,z  are  the  parameters  of  a right-handed  orthogonal  system  of  coordi- 
nates with  +x  in  the  downwind  direction  z to  the  zenith  with  the  origin 
at  the  earth's  surface, 

x,y,T  are  mean  coordinates  of  the  center  point  of  a diffusing  cloud 
t is  time 

a is  a constant  (about  1.25) 
b is  a constant  (about  0.4) 

h source  or  emission  height  (for  Praire  Grass  data  = 46  cm) 
k is  Von  Karman's  constant  (about  0.4) 
r is  a pararieter  of  the  exponential  power  law 
C is  mass  concentration  (mass  per  unit  volume) 

L is  the  Monin-Obukhov  length  scale  (see  appendix) 

Q is  the  source  strength  in  units  of  mass  per  unit  time 
IT  is  the  mean  windspeed 
u*  is  the  friction  speed 

u',  v',  w',  T'  are  instantaneous  wind  and  temperature  fluctuations 

Zo  is  the  surface  roughness  height  (for  Project  Praire  Grass  about 
0,7  cm) 

a.  e»  n.  P»  0),  r,  dummy  variables  or  arbitrary  symbols 
a is  a constant  (about  15) 

C * z/L 

Co  ~ Zo/L 
C = Z/L 
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f,  (|>,  are  universal  functions  of  the  Similarity  Theory 
S is  the  normalized  wind  shear  function 


INTRODUCTION 


This  report  describes  the  concepts  of  similarity  theory  applied  to  con- 
tinuous source  diffusion  problems  in  the  atmosphere. 

Various  papers  have  appeared  expounding  the  use  of  similarity  theory  in 
diffusion  methodology  and  estimation  procedures.  Gifford  [1]  described 
the  similarity  model  for  concentration  at  a downwind  point  due  to  a con- 
tinuously emitting  point  source  up  to,  but  not  including,  a proportionality 
factor.  His  study  indicated  that  the  calculated  power  law  index,  m,  of 
the  inverse  power  law 

Concentration  '(downwind  distance)"'” 

was  in  excellent  agreement  with  data. 

Klug  [2]  developed  empirical  functions  which  relate  similarity  theory  to 
concentration  data  and  found  excellent  agreement. 

Pasquill  [3]  outlined  several  approaches  used  by  others  [1,4,5,]  to  de- 
scribe the  Lagrangian  similarity  approach  to  diffusion.  His  analysis  in- 
dicates that  the  conventional  similarity  hypothesis  fails  to  account  for 
the  effects  of  thermal  stratification  of  the  constant  shearing  stress 
layer  in  the  atmosphere. 

This  report  develops  a proportionality  in  the  concentration  relation  simi- 
lar to  that  reported  by  Klug  [2].  A vertical  distribution  function  is 
selected  for  the  falloff  of  concentration  with  height.  By  use  of  the 
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vertical  distribution,  the  vertical  variance  parameter  is  related  to  a 
similarity  term,  and  estimates  of  cloud  spread  are  given. 

THEORY 

Similarity  Hypothesis  for  Turbulent  Flow 
The  similarity  hypothesis  for  turbulent  flow  in  the  constant  stress  layer 
or  surface  boundary  layer  of  the  atmosphere  is  that  all  statistical  prop- 
erties of  the  flow  depend  only  on  the  ratio  z/L  [6].  For  diffusion,  che 
hypothesis  is  that  the  average  vertical  displacement,  7,  of  a cloud  of 
passive  material  diffusing  throughout  the  surface  layer  is  uniquely  de- 
termined by  the  friction  velocity,  u^^,  and  a universal  function  4i(7/L) 
[1.2,3]. 

This  hypothesis  is  represented  by  the  equation 

^ = bu*(|)(I)  (1) 

where  c = z/L. 


Assuming  that  the  rate  of  increase  of  mean  displacement,  x,  of  the  par- 
ticles of  a diffusing  cloud  is  equal  to  the  average  windspeed,  u,  at  the 
level  z [3], 


dx  _ 
dt 


7(7). 


(2) 


The  expression  for  downwind  maximum  concentration  from  a continuous 
point  source  located  at  the  surface  is  given  by  the  proportionality  [1], 


^axial 


ac 


^u(z) 


(3) 
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Expressions  (1),  (2),  and  (3)  are  the  basic  equations  used  in  the  appli- 
cations of  similarity  theory  to  atmospheric  diffusion  in  the  surface 
boundary  for  a continuous  ground  level  point  source.  These  basic  equa- 
tions are  discussed  below. 

According  to  the  similarity  theory,  u of  Eq.  (2)  is  completely  described 
by  another  universal  function  f(c),  u*,  Zo,  and  Von  Karman's  constant  k. 
As  described  by  [6]: 

u(z)  = -p  [f(c)  - f{Co)].  ^ = z/L,  Co  = Zo/L  (4) 

Therefore,  at  the  level  z,  Eq.  (2)  is  written 


^ = u(F)  = ^ [f(?)  - f(co)]  (5) 

Using  Eqs.  (1)  and  (2)  we  have: 

dx  _ ir(F) 
dz  bu*McT’ 

and  substituting  Eq.  (5)  for  IT  and  setting  T = bkx/L,  ? = F/L,  the 
above  differential  equation  in  dimensionless  variables  F,  F is 


dj  ■ f(F)  - f(co) 


(6) 


When  Eq.  (6)  is  integrated,  the  concentration  function  Eq.  (3)  is 
evaluated  for  a downwind  travel  distance  x where  x = x . 
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There  are  variations  to  the  approach  presented  above.  One  such  variation 
described  by  Panofsky  and  Prasad  [5]  is 


I 


dz  _ b 
dt  a °w 


where  a is  a constant  about  1.25.  Combining  this  with  Eq.  (2)  gives 


^ £ u_  _ a^  _i 

dl  ■ b - b % 


(7) 


Empirical  data  on  o as  a function  of  z and  Zo  can  be  used  with  numerical 
integration  of  Eq.  (7)  to  give  x(z)  [5]. 


Many  papers  on  the  similarity  theory  derive  various  forms  for  the  simi- 
larity functions  <{i  and  f of  Eqs.  (1),  (4).  Tables  A-1  and  A-2  (appendix) 
present  several  of  the  common  functions  for  f and  <p  (the  function  f is 
usually  referred  to  as  the  integrated  wind  profile  function). 


The  method  developed  in  this  report  for  tt.?  relation  (x,  z)  is  yet 
another  application  of  the  similarity  equations  thus  described.  This 
method  consists  of  eliminating  the  need  to  specify  both  f and  <t>. 
Instead,  a second  order  differential  equation  in  (^,  T)  is  derived. 
Let 


S((;) 


k ^ du 


(8) 
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Hence,  f(^)  = j dn.  Furthermore,  let  $ and  S be  related  in  the 

G o 

following  sense; 

(flc)  = S(c)‘^  (9) 

Equation  (6)  then  becomes 

to 

and  the  second  derivative  of  this,  with  the  appropriate  substitutions, 
gives  the  second  order  differential  equation  in  (c.  t),  with  only  the 
function  S(c)  to  be  specified, 

^(4  ins)-  0 (10) 

dt  Vdt  / c 

The  7(7)  function  is  now  the  solution  to  Eq.  (10)  given  S(;)  and  the 
initial  conditions  d^/dc  [to  = 0 and  ^(to)  = 0.  Equation  (10)  can  be 
solved  numerically,  or  for  certain  functions  S(c),  Eq.  (10)  admi ts  of 
elementary  closed  fonn  representations. 

Table  A-3  gives  various  solutions  to  Eq.  (10)  in  closed  form  given  the 
corresponding  S(i;)  function.  Form  A of  Table  A-3  is  the  method  used  by 
Klug  [2],  and  form  B is  one  of  the  methods  described  in  [4].  Both  of 
these  uses  support  the  argument  that  ^ is  given  by  Eq.  (9). 

Particular  f,  iji,  or  S must  be  chosen  on  the  basis  of  the  validity  of 
these  functions  in  use  under  a given  thermal  stratification.  For  example, 
if  the  atmosphere  is  unstable,  L < 0,  a valid  S would  be  the  KEYPS 


8 


relation  or  the  formulation  credited  to  Dyer  [7],  If  the  atmosphere 
is  stable,  L > 0,  then  S could  be  selected  as  the  KEYPS  function,  or 
the  relation  as  in  B of  Table  A-3.  In  the  special  case  where  |Lj  -►  <», 
qeneral  aqreement  is  for  S 1. 


The  basic  startinq  point  of  similarity  and  its  application  to  the  diffu- 
sion problem  is  with  Eqs.  (1),  (2),  and  (3).  Equations  (6)  and  (7)  are 
two  wcivs  tn  describe  the  x(t)  relation.  However,  Eq.  (10)  usinq  the 
identity  Eq.  (9)  will  he  used  in  this  report. 


Replacinq  Proportionality  with  Equality 

Kluq  [2]  presents  in  his  paper  a function  fg(c)  which  describes  the  con- 
centration function  oiven  by 


^axial  ^9^'^ 


(11) 

^ Z U ( 2 ) 

His  analysis  leads  to  the  following  step: 

fq(”)  = r,q(-)  (12) 

wh'Me  at  ^ 0,  corresponding  to  111  - ■ for  neutral  stability,  q(0)  = 1. 


0 I an  then  he  evaluated  according  to 


,,  ^axial  — 
fi  = • z u(z) 


neutral  stability 


given  a set  of  data  on  1 he  appropriate  stability  functions 


for  near  neutral  conditions.  Klun  uses  the  te^t  data  tor  near  neutral 
conditions  from  the  Project  Prairie  liras'-  diffusion  expo?-'  ents  [8].  His 


1 


J 

I 

I 

mean  values  of  G are  as  follows:  = 0.09,  = 0.12,  G2qq|^i  = 

G.nn-  = Gionn„  = 0.13.  The  subscripts  refer  to  the  distances  to  the  vari-  i 

400m  oOOm 

ous  arcs  used  in  the  field  test  program.  Klug  points  out  that  this  be- 
havior of  G can  be  ascribed  to  the  fact  that  the  source  height  in  the 
Prairie  Grass  trials  was  not  at  the  surface,  but  at  sampler  height  z = 

46  cm  above  the  ground  for  most  runs,  and  at  150  cm  above  the  ground  for 
runs  63-68. 
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After  G has  been  determined,  Klug  assumes  that  the  sanu?  set  of  G values 
applies  in  diabatic  conditions,  and  the  function  g(rj  can  therefore  be 
evaUated.  .Jhen  he  computed  the  values  of  g,  there  was  no  functional 
dependence  g(',).  Klug  did  find  a significant  correlation  with  ^/u*. 

His  linear  regression  curve  is  g( ’^/u*)  = -1.28  * l.AJl. 

Klug's  resulting  formula  for  estimating  the  maximum  downwind  concentration 
for  a continuous  point  source  located  at  ground  level  is 

, (-3.2.1  V“‘  " '•‘'8) 

« rro"  ■■■ 

The  use  of  this  estimator  is  reported  to  be  in  excellent  aoreement  with 
the  date  [2]. 

An  fg  function  will  be  developed  in  a similar  manner  in  this  report. 

The  function  fg  relating  observed  concentration  to  the  similarity  ex- 
pression of  tq.  (3)  will  be  expressed  as 

10 

i 

if 


r 


^9(c)  = G(c)g(;) 


(13) 


The  functions  G(c)  and  g(c)  will  be  discussed  below. 

Consider  how  the  effects  of  h,  a nonzero  source  release  height,  is  taken 
into  account.  Sutton  [9]  used  the  method  of  images  to  account  for  the 
effects  of  an  impervious  earth  and  a nonzero  source  height.  Consider 
the  extension  of  his  method  to  a generalized  distribution  of  particles 
in  the  vertical  using  the  functional  form 


x(z)  = exp 


.1^ 


ii: . exp  - JifM: 


(14) 


The  parameter  r is  assumed  constant  together  with  a prescribed  height  of 
emission  h.  B is  a function  of  downwind  distance  x.  The  vertical 
spread  of  particles  is  then  determined  by  the  ratio  x(z)/x(z  ) where 
is  the  height  at  which  the  maximum  concentration  will  occur.  A 

iTla  A 

concentration  function  for  the  xz  plane  will  then  be  constructed  as 


is  a one  dimensional  concentration 

function.  Assume  for  a surface  based  point  source.  Then 

^ zu(T) 


:'x,z)  = 


x("max^ 


g(c) 


1 


Z^u(z) 


. Now  suppose  measurements  are  made  of  con- 


centration using  samplers  located  on  a grid  with  a fixed  sampler  height, 
Zj.  The  comparison  of  the  data  to  the  above  theoretical  form  should  use 


I 


11 


1 


(15) 


C(x,z^) 


x(zj 


' 9(1)  ^ 


ru(z) 


Comparing  Eqs.  (11),  (13),  and  (15)  yields 

x(Zc) 


G(T)  = 


^^^max^ 


(16) 


To  find  set  the  derivative  of  Eq.  (14)  with  respect  to  z to  zero 
max  TV/  r 

and  solve  for  z = For  0 _<  z < “,  the  root  to  = 0 is  an 

absolute  maximum.  With  n = the  solution  for  z is  derived  from 

n max 

f* 

TANHCud  + n)*"  - u(l  - n)"^)  = n,  where  u = 2Cr'~  T)B 
and  TANH  is  the  Hyperbolic  tangent  function. 

Now  B in  Eq.  (17)  can  be  found  by  invoking  the  definition  of  z[l,2,3]; 
that  is,  z is  the  first  moment  of  particle  distribution  in  the  vertical 


— o 

Z = — 


/"  z C(x,z)dz  zx(z)dz 


C(x,z)dz 


J x(z)dz 


(18) 


Using  Eq.  (14), 


z - h6(l/r,  1^)  = [1  -6(2/r,  (19) 

where  r(o)  = /“  n“”'e''^dn  is  the  ganma  function,  and  6(a,e)  * 


0 then 


Note  that  g—  -*■ 


b' '’•(?) 


(20) 


Remember  YCx)  is  known  and  for  x = x B is  a function  of  downwind  travel 
distance,  B = B(x)  (actually  B is  the  vertical  size  parameter  for  the 
diffusing  cloud,  for  if  h = 0,  r = 2,  then  B = 2o^  where  is  a commonly 
used  descriptor  of  vertical  cloud  size  to  be  defined  later).  Inserting 
B(x)  into  Eq.  (17)  defines  z„.^  which  indicates  G((;)  = G(x)  which  was 
shown  empirically  by  Klug  (see  page  10  of  this  report). 


Similarity  Theory  and  the  Vertical  Distribution  z^ 

The  statistic  z^  of  the  distribution  of  diffusing  particles  from  the 
position  they  would  otherwise  occupy  as  a result  of  average  displacement 
7 is  defined  as  follows: 

ao 

F"  = y z^C(x,z)dz //  C(x,z)dz  (21) 

substituting  Eq.  (14)  in  the  C(x,z)  function  as  shown  on  page  11,  then 

r = I z2x(z)dz  / / x(z)dz  = (22) 

o ' ° 

For  simplicity  r = 2 gives  the  familiar  form 

B = 2a?  = 2(P  - h?)  (23) 
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where  the  quantity  is  the  familiar  standard  deviation,  (P  - I )^, 
where  I = h,  of  a normally  distributed  (Gaussian)  probability  function,  but 
is  known  as  the  vertical  size  parameter  of  a diffusing  cloud  in  the 
context  of  this  report. 


Substituting  B = 2a^  in  Eq.  (14)  and  with  r = 2, 


x(z)  = exp[-Js(z  - h)2/o2]  + exp[-5s(z  + h)^/a2T  (24) 

The  above  form,  Eq.  (24),  is  widely  used  in  concentration  predictors  as 
the  vertical  distribution  function  [9]. 


Eval uation 

Now  that  the  diffusion  model  has  been  defined,  constants  and/or  functions 
like  b,  r,  and  g{I)  can  be  evaluated.  First,  the  values  of  b and  r will 
be  found  by  using  data  on  as  defined  in  Eq.  (23)  from  Haugen  et  al. 
[10].  The  predictor  is 

a2  = P - h2  = (25) 
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There  are  48  runs  for  v;hicha  at  v = 130  meters  could  be  calculated. 

z 

The  similarity  functions  for  S(c)  are  as  follows  (Table  A-3): 

(1  - 15  z/L)"'^\  z/1.  < 0 
S(0  = (26) 

(1  + 3 z/L)  , z/L  > 0 

zjl  - Ri 

where  6 = — — k with  z =2  meters  and  R.  is  the  Richardson  number 

Zr/L  1 

for  the  2-meter  level.  Von  Karman's  constant  k is  set  at  0.40.  The 
value  for  r was  chosen  to  be  1.5  for  the  following  reasons.  Elliot's 
[11]  analysis  of  Prairie  Grass  data  for  the  vertical  concentration 
distribution  did  not  reveal  any  substantial  variation  of  shape  with 
stability,  and  his  average  r was  1.49.  The  value  of  b in  Eq.  (1)  was 
adjusted  to  obtain  the  best  fit  of  Eq.  (25)  to  Haugen's  data.  A value 
of  b = 0.48  was  obtained. 


The  resulting  value  of  bk  = 0.19  is  not  far  from  the  experimental  value 
of  0.18  obtained  by  Pasquill  from  these  same  data  [3].  Figure  1 shows 
the  comparison  of  theory,  i.e.,  Eq.  (25),  to  the  data.  Very  good  agree- 
ment is  obtained  for  all  L such  that  1/L  < 0.03.  There  is  fair  agreement 
for  the  very  stable  condition  1/L  > 0.03. 


The  function  g(c)  can  now  be  determined  from  the  data  by  the  relation 
9 " \(x]^  where  is  the  measured  field  of  concentration  for 


a run  normalized  by  the  source  strength  for  that  run.  Fifty  runs  were 
obtained  from  approximately  70  tests  of  the  Prairie  Grass  Program  to 
utilize  in  evaluating  g.  As  with  Klug  [2],  no  relation  such  as  g g(c) 
was  found.  Instead,  after  finding  an  average  g,  g,  for  each  run,  50  g 
values  were  correlated  with  o„/u*.  The  linear  regression  is 

gCo^u*)  = -0.035o^/u*  + 0.186  (27) 

This  regression  line  and  the  data  are  plotted  in  Figure  2.  To  compare 
this  linear  regression  line  to  that  of  Klug's,  Eq.  (27)  has  to  be  scaled 

g(o„/u*) 

by  1/0.13.  Ilith  this  scale  — — ""  -0.27  o^/u*  + 1.43.  The  compari- 
son is  in  good  agreement  considering  the  different  approach  taken  by 
the  two  methods.  (See  page  10  this  report.) 


Our  prediction  formula  for  concentration  in  the  xz  plane  from  a con- 
tinuous point  source  for  a near  ground  level  release  is 


, . (-0.035  oju^  + 0.186) 

i 

x^'max^  tu{T) 


(28) 


The  data  of  Tables  A-5,  A-6,  and  A-7  depict  the  applicability  of  the 
similarity  theory  to  the  diffusion  problem.  The  data  on  comparisons 
in  Table  4-5  show  very  good  agreement  with  the  overall  ratio  of  predicted 
(z^  - h^)*/^  to  observed  0^  val ues  equal  to  1.04.  The  concentration  com- 
parisons of  Table  A-6  give  good  agreement.  There  is  an  anomalous  value 
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of  the  800-meter  value  of  concentration  for  run  39.  But  the  overall 
ratios  of  predicted  to  observed  concentrations  are 


50  = 1-02 

r^oo  = 1 -06 

•^800  ■ 

100  ^ 

•"400  =1-21 

The  comparison  of  the  vertical  distribution  Eq.  (14)  using  r = 1.5  is 
given  in  Table  A-7.  Here  Elliot's  [12]  data  on  the  ratio  of  vertical 
concentration  at  height  z to  that  at  height  z = 0.5  is  used.  The  distri- 
bution Eq.  (14)  does  an  excellent  job. 

CONCLUSIONS 

The  Similarity  Theory  aided  with  empirical  derivations  has  been  shown  to 
be  applicable  to  the  continuous  point  source  problem  of  atmospheric  dif- 
fusion. The  stability  functions  f,  4),  and  S allow  the  theory  to  cover 
all  atmospheric  stabilities  of  neutral,  stable,  or  unstable. 

This  report  has  shown  that  the  solution  to  a second-order  ordinary  dif- 
ferential equation  in  (c,  c)*  Eq.  (10),  which  needs  only  a function  for 
S(c),  the  nondimensional  wind  shear, specifies  the  change  of  z versus  x. 
(Note  that  the  integration  can  be  carried  out  even  if  S is  in  the  form  of 
data  points  directly  computed  from  a wind  profile,  although  this  was  not 
done  in  this  report.) 


17 


The  solution  to  Eq.  (10)  in  z and  x is  then  used  to  compute  the  concen- 
tration C.  This  report  has  also  developed  a proportionality  for  Eq.  (3) 
based  on  an  analytical  treatment  of  the  vertical  diffusion  term  (Eq.  14). 

The  vertical  size  parameter  z^  has  been  given  a functional  form  (Eq.  22) 
and  has  been  shown  to  be  in  excellent  agreement  with  data  over  all  atmo- 
spheric stabilities. 

A restraint  on  h,  the  source  emission  height,  is  that  it  should  not  be 
greater  than  a few  meters  above  ground  level.  For  an  excellent  suimary 
of  restrictions  on  h and  x in  similarity  theory  see  Pasquill  [16],  sec- 
tion 3.3. 
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APPENDIX 


CALCULATION  PROCEDURES 

The  calculation  procedures  for  finding  L,  u*,  6(R^),  R^- . and 
merical  calculation  of  S(c)  from  profile  data)  are  as  follows: 


R.  = g/T 


Z AlnZ  where 


T = °c  + 273.16 


r = 0.0098  deg  k/m 


g = 9.8  m/sec 


Zg  = geometric  mean  of  the  profile  layer 
V is  the  wind  profile  data,  m/sec 


In  stable  flow. 


_ kg  dG/aZ  _ kg  AO  i /i  _ ll/l 

■ ITg  aWz"  - u~W 


Tm  = [1  ■* 


AV^  - AV2 
AV^  - 2AV2 


] = 


u*  = k^s"^  AV/AlnZ 


8 = R. 


1 AV^  - AV^ 

ivr~" 


In  unstable  flow. 


Z/L  = R, 


3^^ 

u*  pOd 


i kZgH 


1/L  = 


T R.Z. 


I Z, 


(a  nu- 
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H = CppK 


L© 
H aZ 


ku^Z 

H 


AVif  - AVj 
- 2aVJ 


u*  = k(j>*^  AV/AlnZ 


Y = 15 
where 


Cp  = 2.24  cal/g  deg 

p = 0.0012  g/cm^ 

k = 0.4 

0 = ° 


C + 273.16 


TABLE  A-1 


INTEGRATED  WIND  PROFILE  FUNCTIONS  COMMONLY  IN  USE 

A.  KEYPS  (Yamamoto  [13]) 

f(i;)  = 5-  2TAN'6-2  TANH“^5,  where  6**  - - 1 = 0 with  ^ 

about  7,  c = z/L 

B.  Moni n-Obukhov  [6] 

In;  + Be,  B > 0 

f(0  = 

In;  B = 0 

; - R. 

Bl.  As  above  but  with  b variable  with  b = 5 — • 

;Ki 

; = z^/L,  z^  a reference  height, 

R.  = Richardson  number  at  level  z 
1 r 

C.  Businger,  Dyer  [14] 

f(;)  = In  2 TAN-Ie 

where  t = (1  - y about  15 
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TABLE  A-2 


THE  UNIVERSAL  FUNCTION  *(;)  OF  EQ.  (1) 


A.  Monin  (1959)  [15] 

/ d f ^ 1/4 

♦ (<:)  = [1  - (^1  ] where  f is  the  integrated  wind  profile  function. 

B.  Gifford  [1] 

,U)  - 1/2  t 1).  . > 5 

I T aC 

C.  KOEHLER  (1967)  [16] 

(|i(c)  = (1  + Be)"',  6 > 0 
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TABLE  A- 3 


THE  RELATION  OF  x AND  z USING  EQ.  (10)  AND  THE  FOLLOWING  S(;)  FUNCTIONS 
(using  the  initial  conditions  dc/dcico  = 0,  c(co)  = 0). 

A.  KEYPS  [13] 

S(i;)  given  by  S'*  - a'c  S^  - 1 = 0 
then 

I = . (S'*  - 2j2  - 3)tan'^S  - 1/2(S‘*  - 3)(So  - 2tan"^So) 

(1  + So)(i  - s)-| 

+ 1/2(S‘*  + 2S2  - 3)ln  (1  . Sord'  VsyJ 

where  So  - aco  So  - 1 = 0 

B.  Monin-Obukhov  [4] 

S = 1 + e? 

C = {12£lnc/co  + 6e2lnc/Co  + 12(e  - Eo)  + 3(e  - 3Eo)(e  - Eo) 

+ 2(e  - eo)^{2e  - Eo)}»  where  e = 6^,  eo  = Bto 

C.  Businger,  Dyer  [11] 

S = (1  - 

_ P (1  - e^)(l  + Eo^)  (1  + e)(1  - Eo) 

« = 3?  (1  + Ei!)(1  - Eo^)  ^ (1  - E)rT+"Eo) 

+ E3(tan*^Eo  - tan'^s)} 

-i  -1  " 

where  e » S , e = So 
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TABLE  A-4 


VALUES  OF  G’*(x)  FOR  TYPICAL  STABILITY  CONDITIONS 


(Zo  = 0.007 

m,  b = 

0.48,  k 

= 0.4,  r 

= 1.5,  h = 

0.46  m 

L 

Distance 

From  Source,  x 

Val  ue 

Stability 

Used 

50  m 

100  m 

200  m 

400  m 

800  m 

Extremely  Unstable 

L > -5 

-1 

1.07 

1 .02 

1 .00 

1 .00 

1 .00 

Moderately  Unstable 
-100  1 L < -5 

-50 

1.27 

1 .04 

1.01 

1 .00 

1 .00 

Near  Neutral 
|L|  >100 

-200 

1.30 

1.12 

1.05 

1 .02 

1 .00 

+200 

1 .35 

1.15 

1 .07 

1 .03 

1 .01 

Stable 

10  < L < 100 

+50 

1.42 

1.79 

1.09 

1.05 

1.02 

Very  Stable 

L < 10 

+1 

2.17 

1 .63 

1.36 

1 .26 

1 .13 
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TABLE  A- 5 


COMPARISON  OF  VERTICAL  SPREAD  [z2-h2]V2  at  100  tn. 
a OBSERVED  (HAUGEN,  ET  AL.),  [z2-h2]V2  PREDICTED, 
AND  R,  RATIO  OF  PREDICTED  TO  OBSERVED 


Run 

No. 

o(ni) 

[z2.h2]i/2 

R 

Run 

No. 

o(m) 

[z2-h2]V2 

R 

14 

1.4 

1.18 

.74 

41 

3.4 

3.96 

1.16 

15 

7.7 

7.46 

.97 

42 

4.4 

4.08 

.93 

16 

10.1 

8.86 

.88 

43 

7.0 

6.60 

.94 

17 

3.7 

3.57 

.96 

44 

6.3 

5.79 

.92 

18 

2.7 

3.30 

1.22 

45 

4.5 

5.04 

1.12 

19 

6.9 

6.01 

.87 

46 

3.3 

4.00 

1.21 

20 

5.8 

5.37 

.93 

48 

5.8 

5.2 

.90 

21 

4.2 

4.32 

1 .03 

49 

4.5 

5.69 

1.26 

22 

4.6 

4.32 

.94 

50 

5.9 

5.67 

.96 

24 

4.6 

4.38 

.95 

51 

4.8 

5.54 

1.15 

25 

9.0 

8.65 

.96 

52 

10.0 

7.44 

.74 

26 

5.7 

5.6 

.98 

54 

3.9 

3.58 

.92 

27 

5.8 

5.54 

.96 

55 

4.4 

4.17 

.95 

29 

3.4 

3.46 

1.02 

56 

4.4 

3.91 

.89 

30 

5.6 

5.33 

.95 

57 

4.9 

4.85 

.99 

31 

6.5 

5.07 

.78 

58 

1.5 

2.75 

1.83 

32 

1.6 

2.18 

1.36 

59 

2.0 

2.90 

1.45 

33 

5.1 

5.17 

1.01 

60 

3.9 

3.80 

.97 

34 

5.3 

5.14 

.97 

61 

4.8 

5.60 

1.17 

35s 

3.4 

3.74 

1.10 

62 

6.9 

5.73 

.83 

36 

1.9 

2.54 

1.34 

65 

3.4 

3.28 

.96 

37 

4.5 

4.02 

.89 

66 

2.1 

2.70 

1.29 

38 

3.8 

4.06 

1 .07 

67 

3.3 

3.49 

1.06 

40 

3.4 

2.65 

.78 

68 

1.9 

2.72 

1.43 

average  ratio  = 1.04 
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TABLE  A-6 


VALUES  OF  CALCULATED  PARAMETERS  FOR  ALL  PRAIRIE  GRASS  RUNS, 
AND  RATIOS  OF  PREDICTED  CONCENTRATIONS  BY  THEORY 
TO  MEASURED  CONCENTRATIONS  FOR  THOSE  RUNS  SELECTED  BY  THIS  STUDY 


Run 

No. 

L 

(R^)2m 

*^50 

O 

o 

*^200 

^^400 

*^800 

1 

-14.7 

.197 

★ 

2 

-27.0 

.128 

ir 

3 

0.35 

.020 

.5 

4 

0.35 

.020 

.5 

5 

-32.3 

.434 

* 

.93 

.93 

1.03 

1.31 

1 .46 

6 

-111. 

.429 

★ 

1.08 

1 .15 

1.27 

1 .52 

1.65 

7 

-14.1 

.336 

★ 

8 

-25.6 

.333 

* 

.46 

.47 

.48 

.68 

.91 

9 

-19.6 

.434 

* 

1.04 

1 .05 

1.07 

1.17 

1.32 

10 

-11.0 

.358 

★ 

.69 

.75 

.65 

.75 

11 

-90.9 

.546 

★ 

.95 

.93 

.76 

.76 

.81 

12 

-58.8 

.586 

★ 

1.16 

1.13 

.91 

.85 

.64 

13 

0.35 

.020 

.5 

14 

0.35 

.020 

.5 

15 

-8.13 

.253 

* 

.52 

.53 

.50 

.53 

.88 

16 

-4.59 

.267 

* 

17 

-56.0 

.209 

* 

.84 

.84 

.98 

1.21 

1.46 

18 

31.4 

.213 

* 

.94 

1.06 

1.24 

1.45 

1.69 

19 

-22.2 

.460 

* 

20 

-50.0 

.680 

* 

1.15 

1.30 

1.24 

1.43 

1.65 

21 

242. 

.425 

.007 

.74 

.76 

.85 

.93 

1.02 

22 

238. 

.507 

.008 

.80 

.79 

.80 

.89 

1.05 

23 

259. 

.414 

.007 

.97 

1.00 

1.05 

1.18 

1.28 

24 

302. 

.397 

.006 

1 .09 

1.13 

1.15 

1.14 

1.25 

25 

-5.00 

.222 

* 

26 

-34.5 

.474 

* 

1.22 

1.18 

1.09 

1.28 

1.64 

27 

-37.0 

.459 

* 

1.20 

1.29 

1.14 

1.07 

1.01 

28 

26.8 

.167 

.043 

1.24 

1.38 

1.84 

2.15 

2.17 

29 

40.4 

.273 

.033 

30 

-52.6 

.504 

* 

31 

-100. 

.572 

* 

1.18 

1.34 

1.43 

1.67 

1.24 

32 

6.31 

.094 

.102 

1 .81 

1.17 

1.29 

1.40 

1.40 

33 

-76.9 

.547 

* 

.92 

.87 

.83 

1.01 

1.36 

34 

-83.3 

.679 

* 

1.12 

1.03 

.92 

.89 

.82 

35 

2.20 

.051 

.224 

35s 

65.5 

.253 

.023 

.88 

.96 

1.15 

1.38 

1.53 

36 

9.21 

.084 

.091 

1.61 

1.29 

1.68 

2.59 

1.72 

37 

114. 

.315 

.014 

1.03 

1.16 

1.45 

1.86 

2.94 

38 

129. 

.295 

.013 

.82 

.74 

.77 

.80 

.87 

30 


Table  A-6 (Cont) 


Kun 

No. 

L 

U*  (R^)2m 

«50 

o 

o 

^200 

o 

o 

*^800 

39 

7.32 

.102  .104 

1.33 

1.92 

2.59 

3.27 

40 

10.0 

.101  .086 

41 

105. 

.320  .015 

.60 

.55 

.60 

.57 

.56 

42 

136. 

.421  .012 

.93 

.99 

1.18 

1.75 

2.02 

43 

-13.5 

.395  * 

.89 

.99 

.89 

1.09 

.98 

44 

-27.8 

.443  * 

1.02 

1.08 

1.08 

.96 

1.16 

45 

-125. 

.430  * 

1.06 

.99 

.90 

.91 

.77 

46 

112. 

.372  .015 

.85 

1.06 

1.40 

1.85 

2.89 

47 

-4.67 

.269  * 

48 

-7.14 

.564  * 

1.30 

1.21 

1.21 

1.13 

1.04 

48s 

-5.55 

.257  * 

49 

-31.3 

.494  * 

1.19 

.98 

.98 

1.20 

1.19 

50 

-32.3 

.493  * 

.96 

1.02 

1.29 

1.41 

2.29 

51 

-37.0 

.498  * 

1.04 

1.09 

1.20 

1.54 

2.93 

52 

-8.33 

.344  * 

53 

5.37 

.103  .129 

1.45 

1.26 

1.43 

1 .83 

2.10 

54 

47.1 

.271  .030 

.84 

.85 

.96 

.95 

.99 

55 

166. 

.418  .010 

.90 

.87 

.94 

1.19 

1.28 

56 

91.0 

.330  .018 

.85 

.90 

1.04 

1.11 

1.14 

57 

-333. 

.525  * 

1.11 

1.33 

1.11 

1.22 

1.31 

58 

18.0 

.109  .053 

1 .02 

.65 

.61 

.63 

.64 

59 

13.3 

.150  .075 

.80 

.56 

.59 

.55 

.62 

60 

70.6 

.320  .022 

.78 

.77 

.77 

.93 

.93 

61 

-34.5 

.572  * 

1.33 

1.25 

1.07 

1.07 

1.23 

62 

-29.4 

.381  * 

.96 

.84 

.69 

.62 

.63 

63 

0.35 

.020  .5 

64 

0.35 

.020  .5 

65 

59.1 

.299  .025 

1.20 

.89 

.87 

.87 

1.02 

66 

23.9 

.188  .049 

67 

84.7 

.299  .019 

1.14 

.85 

.80 

.84 

.92 

68 

24.3 

.160  .048 

1.05 

.92 

1.04 

1.22 

.97 

overall 

average  ratios 

1.02 

1.00 

1.06 

1.21 

1.29 

♦For  unstable  conditions  z/L  * 


31 


TABLE  A- 7 


COMPARISON  OF  RATIOS  OF  GAS  CONCENTRATIONS  AT  HEIGHT  Z 
TO  CONCENTRATION  AT  O.S  *'i.  OBSERVED  Ro . AND  PREDICTED  Rp 


AT 

100  M FROM 

RELEASE 

POINT 

(R„ 

VALUES  FROM 

ELLIOT 

[5]) 

Run 

Z 

Run 

Z 

Run 

Z 

No. 

cm 

R. 

R 

P 

No. 

cm 

Ro 

R 

P 

No. 

cm 

R„ 

R 

P 

17 

100 

.95 

.93 

18 

100 

.93 

.92 

19 

100 

.96 

.97 

150 

.87 

.85 

150 

.82 

.83 

150 

M 

.93 

250 

. 70 

.67 

250 

.61 

.63 

250 

.73 

.83 

450 

.42 

.35 

450 

.25 

.30 

450 

.52 

.61 

750 

.14 

.096 

750 

.022  .068 

750 

.30 

.34 

.’0 

100 

.98 

.96 

21 

100 

.94 

.95 

22 

100 

.94 

.95 

150 

.96 

.91 

150 

.87 

.88 

150 

.88 

.88 

250 
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